Interplay between charge-lattice interaction and strong electron correlations in 
cuprates: phonon anomaly and spectral kinks 
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We investigate the interplay between strong electron correlations and charge-lattice interaction 
in cuprates. The coupling between half breathing bond stretching phonons and doped holes in the 
t-t'-J model is studied by limited phonon basis exact diagonalization method. Nonadiabatic electron- 
phonon interaction leads to the splitting of the phonon spectral function at half-way to the zone 
boundary at q 3 = {(±7r/2, 0), (0, ±7r/2)} and to low energy kink feature in the electron dispersion, in 
agreement with experimental observations. Another kink due to strong electron correlation effects 
is observed at higher energy, depending on the strength of the charge-lattice coupling. 
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There is a growing confidence that strong electron- 
phonon interaction (EPI) manifests itself both in vibra- 
tional [l[ and electronic Q spectra of cuprates. The most 
puzzling feature of the cuprate phonon spectra is the 
anomaly of the half breathing bond stretching (HBBS) 
phonon occurring at half-way to the Brillouin zone (BZ) 
boundary in the [100]-direction, while the most debating 
feature of the electronic spectra are the kinks observed 
in Angle Resolved Photoemission Spectra (ARPES). 

It was recently realized in the experimental commu- 
nity that it is highly important to measure both phonons 
spectra and ARPES on the same sample just to verify 
possible link between HBBS phonons and ARPES |3(. 
The results of the above studies support close connec- 
tion between HBBS phonon anomaly and lowest energy 
kink in ARPES 0] and, thus, a model describing both 
anomalies within the same approach is strongly needed. 

In the present Letter we study the low density limit 
(one hole on 4 x 4 lattice) of the extended t-t'-J model 
where holes are coupled to HBBS phonons Q. In or- 
der to calculate the phonon spectral function (PSF) and 
hole spectral function (HSF), we generalize a recently 
introduced approach [f|, based on limited phonon basis 
exact diagonalization (LPBED), without adopting self- 
consistent Born (SCBA) and spin wave approximations 
@, 0] ■ This method treats the non-adiabatic effects of 
the quantum phonon very effectively without approxi- 
mations for the magnetic degrees of freedom. Because of 
exponential growth of the basis with size of the system, 
the 4x4 lattice has considerably denser quantum states 
than the smaller \/T0 x \/To system Q , so that it is pos- 
sible to resolve fine structure of the PSF and HSF. For 
the first time we are able to reveal the shape of the PSF 
while previous studies were restricted to at most 

the second moment of the response. 

We show that EPI can lead to the splitting of the PSF 
at half-way to the BZ boundary in the [100] direction. We 



demonstrate that the splitting can be easily smeared out 
by a very small broadening and, thus, the mysterious eva- 
sive behavior of double-peak structure can be attributed 
to tiny variations of chemical composition, crystal quality 
and/or experimental setup. We argue that the splitting 
is a rather general phenomenon arising when the phonon 
branch interacts with a rather soft electronic excitation 
and we show that the HBBS phonon anomaly is linked to 
the lowest energy kink observed in the ARPES. Finally 
we emphasize that the same model supports the spec- 
tral kink at higher energy referred to colloquially as the 
waterfall 
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McQueeney et al. [13( observed anomalous lineshape of 
HBBS phonons around q = (n/2,0, 0) and reported its 
strong temperature dependence. Subsequent studies con- 
firmed that the HBBS phonon anomaly is due to the cou- 
pling of the HBBS phonons to the doped carriers. Indeed, 
the anomaly is absent in undoped compounds til and the 



hardening of the phonon spectra with heating [13l Il4| ex- 
cludes such sources of anomaly as anharmonicity or struc- 
tural inhomogeneity. It is also agreed that HBBS phonon 
anomaly is located around q = (tt/2, 0, 0) at any doping. 
On the other hand, there is a controversy on the rela- 
tion between the HBBS phonon anomaly and the lowest 
energy kink in the ARPES. This relation is often denied 
although the measurements of the phonon spectra and 
ARPES on the same sample of Bi2Sri.6Lao.4Cu206+a; 
suggest that the softening of the HBBS phonon mode 
matches the energy and momentum of this kink Q . An- 
other much debated question is the interpretation of the 
PSF structure at HBBS phonon anomaly. Actually the 
phonon peaks at q s are poorly defined and two scenar- 
ios are possible. According to the single-branch inter- 
pretation (SBI) there is softening and a very large in- 
trinsic linewidth of a single peak near q = (tt/2, 0,0). 
On the contrary, the double-branch interpretation (DBI) 
implies the presence of two peaks: the first one is lo- 
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cated at the unrenormalized energy and the second one 
is broad and soft. The pioneering study in Ref. [l~3| was 
followed by [l7| whose results supported SBI. Next stud- 
ies of similar material Lai.g75Bao.i25Cu04 led authors to 
conclusion about relevance of DBI [l5|. Likewise, mea- 
surement of HBBS in YBa 2 Cu306+a; was interpreted in 
the framework of DBI [l~8[. Recently a "contamination 
hypothesis" [l|| has been suggested, claiming that the 
SBI is correct because the inelastic neutron scattering 
measurement picks up the intensity from the k-vicinity 
(ir/2, k,0) of q = (ir/2,0, 0) point and, thus, the "nor- 
mal" component comes from transverse phonons. How- 
ever, the above reasoning contradicts the results obtained 
on La 2 - x Ba 2; Cu04 + «5 at i = 0.14 ± 0.01 by in- 
clastic x-ray scattering technique which has higher mo- 
mentum resolution than that in neutron scattering mea- 
surements. " Contamination hypothesis" implies that the 
"normal" component must disappear by impro ving mo- 
mentum resolution. Instead, the results of Ref. [20j sup- 
port DBI. One can guess that the SBI vs DBI depends 
on the chemical composition Ba-doped 15, H(j vs 
Sr-doped 14|,ll7| compound. Besides, since the softening 
is very sensitive to the broadening caused by decay chan- 
nels and/or experimental resolution, SBI vs DBI must 
be strongly dependent on the compound, sample quality, 
and experimental setup. 

We will show that many of the questions discussed 
above can be understood within the t-t'-J model includ- 
ing the coupling with HBBS. On the contrary density 
functional theory calculations [21 , 22 do not predict any 
HBBS phonon anomaly [23|. Other approaches, espe- 
cially those associating the anomaly with stripes, have 
difficulty with the position of HBBS phonon anomaly. 
Scenario suggested in Ref. [24j relates the HBBS soften- 
ing to the Kohn anomaly at double Fermi momentum 
2fcp along the Fermi surface of stripes 24|. Here, in con- 
trast to experiment, the softening must be 9- independent. 
In a different theoretical proposal, HBBS anomaly is as- 
sociated with stripe mediated collective charge excita- 
tions [25I ] or incommensurate low energy spin-fluctuations 
[jjj . Within these scenarios, in contrast with experi- 
ment 14. 15. 2ol 27 , 28| , wave vector of phonon anomaly 
strongly depends on doping level. Finally we note that 
phonon softening [Tol. broadening Q , and correct po- 
sition of the anomaly [11| have already been qualitatively 
explained by the coupling between phonons and density 
response of the t-J model although none of those studies 
considered the shape of the phonon spectral function. 

The Hamiltonian of the t-t'- J-Holstein model in 2D is 
the sum of the electronic part and hole-phonon coupling 
Hamiltonian 

Htt'j = ~ * C l+S,a c i,<? ~ * C j+S',a Ci >< 7 
i,S,a i,8' ,o 

+ ^Y S ' l+sSl - , (1) 



H' = uj a \.n a w + (Mq^e^ 1 ** (1 - ni)a qill + H.c. 

q,li i,q,n 

Here J is the exchange interaction constant of the spin- 
spin interaction, t and t' arc hopping amplitudes to near- 
est and next nearest neighbors. At site i (Cu atoms), 
Si is the i-spin operator, Ci iCr is the fermionic operators 
in the space without double occupancy, and n, is the 
number operator. a q ^ is the phonon annihilation oper- 
ator with momentum q and fx = x or y indicates the 
longitudinal polarization of the oxygen vibrations along 
the direction of the nearest neighbor Cu atoms (Oxy- 
gen atoms are located at i?,- + a/2/i, a being the lat- 
tice parameter). M q ^ is the matrix element of the EPI: 
M q fl = qijJqI \f~N2i sin((j^/2) where N is the number of 
lattice sites and uiq is the frequency of dispcrsionless op- 
tical phonon. The strength of EPI is characterized by 
dimensionless coupling constant A = J2qp, |Mj.^| 2 /4wot. 
We chose parameters which correspond to hole doped 
cupratcs: J = OAt, t' = — 0.25i, ljq = 0.15t, and set 
t = a = H = 1. The PSF is expressed as 
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D(q, lj + irj) = 3 — — j — - — 

7r D (q,ui + irj) — + irj) 

where the self-energy E(<f,u; + irj) is given by 



Z(q,u} + ir)) 
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Here Do (<f, u + irj) is the bare phonon Green function and 
n(g, uS) = P(q, to + irj) + P(q, -co — irj) is the polarization 
insertion with 



P(q,cj + iri) = (^ |O+- 
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We choose the ground state (GS) \ipo) as a linear super- 
position with equal weights of the 4 degenerate states 
corresponding to k = (±^,±1) with energy Eq, rj is 
a broadening factor that shifts the poles of D{q,ui) in 
the complex plane, and O = '^2 i e lq ' Ri (l — rii). The in- 
crease of broadening factor 77 has the physical meaning 
of phonon damping or limited experimental resolution. 

The ground state \ipGs) an d the function P(q,ui + irj) 
are obtained by modified (29j and standard Lanczos 
methods, respectively, within the LPBED method Q in- 
troduced for the t-t'- J-Holstein model. 

In Fig. 1 we plot the PSF D(q, u>) for different wavevec- 
tors along the (1,0) direction in the BZ. For the cho- 
sen values of the model parameters the system under- 
goes, in agreement with [301 ] . a crossover towards strong 
EPI regime for A c ~ 0.5. The anomalous softening of 
the phonon mode at q s = (n/2,0) is already observed 
at moderate values A < A c of the hole-phonon coupling 
(A = 0.3). We stress that for a large value of broadening 
factor T) = 0.08 the phonon peak softens and broadens at 
q s (Fig. la), supporting, thus, SBI. On the other hand, 
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FIG. 1: (Color online) (a) The PSF for A = 0.3 at three dif- 
ferent values of q (q — (0, 0) (dashed green line), q = (n/2, 0) 
(solid red line), and q — (n, 0) (dotted blue line)) with 
r\ — 0.08; (b) the PSF with a smaller broadening factor 77, 
?7 = 0.03. The vertical lines indicate the bare phonon fre- 
quency UJo ■ 
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FIG. 2: (Color online) HSF A(k,u), for A = 0.3 at two 
different values of k (k = (n/2,n/2) (dashed green line), 
k = (0,7r/2) (solid blue line)) with 77 = 0.03. The energy is 
measured with respect to the GS energy. The phonon Green 
function at q = (ir/2, 0) (with the same value of rj) is plotted 
too for comparison (red dotted line). In the inset sketch of 
one of transitions involved in the PSF calculation. 



reducing the broadening factor to 77 = 0.03 two-peak 
structure turns out(Fig. lb), that is in agreement with 
the experimental observations reported in Ref. [Til, [2(| . 

To clarify the physical nature of the splitting, we study 
the hole spectral weight function A(k, to) at the wavevec- 
tors kf that are reached starting from GS wavcnumber at 
(±7r/2, ±7r/2) through phonon momentum q s = (jr/2,0). 
For example, as shown in Fig. 2, by starting from k = 
(— 7r/2,7r/2) we get kf — (0, 7r/2) after absorbing (emit- 
ting) a phonon with momentum — q s (q s ). The lowest en- 
ergy peak in D(q s ,ui) (dotted red curve in Fig. 2) is close 
to the lowest energy peak of A(kf,ui) (solid blue curve 
in Fig. 2). Measuring energy from the GS, denoted by 
vertical line in Fig. 2, one can see that the lowest peaks 
in phonon and hole spectral functions are softer than the 
phonon energy wo=0.15t. We note that the significant 
renormalization of the hole spectral weight, with respect 
to GS one, indicates strong coupling between phonon at 
q s and hole at kf. The high energy peak in D{q s ,u>) is 



located at energy close to luq, and no peak at the same 
energy is found in A{k, uS). It is shown below that it is 
due to the phononic nature of the high energy resonance 
of the PSF. 

To give a simple explanation of the above scenario 
one has to realize that the excited electronic state 
I 7 / , )(ott/2) |0) with the momentum (0,7r/2) and without 
phonons is linked by the matrix elements of EPI with a 
group of 8 degenerate states \ip)°(± 7T /2 ±7r/2) a \.^ |0) where 
the electronic subsystem is in the GS and one phonon 
is excited. The momentum conservation (±7r/2, ±7t/2) + 
q = (0,7r/2) determines the phonon momenta q. The 
energies of the electronic subsystem at corresponding 
momenta are e°(0, 7r/2) and e°(7r/2, 7r/2), respectively. 
Then, choosing a set of the model parameters appro- 
priate for cuprates, the resonance relation e°(0,7r/2) « 
e°(w/2, 7r/2) + ojo is satisfied. Hence, even small matrix 
elements induce strong effects both in electronic and vi- 
brational subsystems. Analytic diagonalization of this 
degenerate 9x9 matrix gives 9 levels. The lowest state 
L has energy below luq and has components on both one- 
phonon and zero-phonon states. This state corresponds 
to the peaks of electronic A(k = (0,7r/2),w) and bosonic 
D(q = (7r/2, 0),w) spectral functions at u < 0.1 (see the 
area around letter L in Fig. 2). In the highest energy 
range, (see the area around letter H in in Fig. 2), both 
spectral functions collect contributions at uj > 0.2. How- 
ever, the iJ-state is of predominantly electronic origin 
with a large peak in the hole spectral function while it 
generates a weak structure in the PSF, hardly observ- 
able in the numerical data. The last 7 degenerate levels 
with energy e°(7r/2, 7r/2) + loq have no projection on the 
vacuum boson state: this explains why the hole spectral 
weight function does not show any peak at this energy, at 
a distance cjo from GS energy. The energies of all levels 
of the above simple analytic solution are in qualitative 
agreement with that provided by the numeric LPBED 
method in Fig. 2. So, we conclude that the doubling of 
the phonon peak is due to coupling between holes and lat- 
tice, that lifts a degeneracy and produces one additional 
state with energy less than loq. 

Generically, the EPI involving electronic states with 
energies considerably larger than the phonon frequency 
does not lead to any splitting since the adiabatic ap- 
proximation is valid and electron produces only a renor- 
malization of the adiabatic potential [3l| . On the other 
hand, when the electronic excitation is soft and its energy 
is comparable with phonon frequency, the nonadiabatic 
corrections play a crucial role [32| and one can observe 
exotic spectral functions as it is seen in Fig. lb and in 
experiments on cuprates. 

The clear connection, established by the above consid- 
erations, between phonon at q s and hole at kf encouraged 
us to search for an interplay between the HBBS phonon 
anomaly and the lowest energy kink in the dispersion of 
quasiparticles in cuprates (33j . In Fig. 3a we plot hole dis- 
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electron correlations and charge lattice interaction are 
crucial in understanding cupratc experimental features. 
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FIG. 3: (Color online) (a) Energy vs (k — kp)/n dispersion 
relations along nodal direction from kF to V point at A = 0.3 
(solid) and A = (dashed). Arrows mark kink positions, (b) 
Dispersion relation as a function of oxygen isotope exchange 
( 16 Q (solid) and 18 Q (dashed)) at A = 0.3. 



persion derived from momentum distribution curves, i.e. 
we find the k for which the hole Green function A(k,u>) 



has maximum at fixed 
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The energy, measured 
from GS, is plotted versus the wavenumber I k — kp 



where kp is the momentum corresponding to the mini- 
mum of the hole dispersion relation. We find that, when 
Ek — Ek F is about wo = 0.15, the curve exhibits a slope 
change related to the coupling between the bare phonon 
energy and the electronic band in the bare t-t'-J model. 
This kink does not appear in absence of the hole-phonon 
coupling. On the other hand, we obtain another kink at 
higher energy (see Fig. 3a) [12}. This distinctive feature 
is related to the strong electron correlations since it is 
observed also in the bare t-t'-J model. At A = the kink 
is located around the exchange interaction e nerg y J, in 
agreement with results by Chakraborty et al. [35j within 
the Hubbard model. Furthermore, the kink is shifted 
at higher energy by increasing the hole-phonon coupling. 
We also investigated the effect of the oxygen isotope sub- 
stitution Q on ARPES. The results (Fig. 3b) show a neg- 
ligible (below O.Oli) shift upon oxygen isotope exchange, 
that is in agreement with recent experimental observa- 
tions by Douglas et al. [36|. All these data point out 
that the scenario based on the interplay between strong 
electron correlations and hole-phonon interaction is able 
to capture many physical distinctive features of high tem- 
perature superconductors. 

In conclusion, we showed that the EPI, in the presence 
of strong correlations, can lead to the splitting of the 
phonon spectral function at half-way to the BZ bound- 
ary in the [100] direction. We demonstrated that the 
splitting can be easily smeared out by very small broad- 
ening of the eigenstates. The same physical mechanism 
can explain both HBBS phonon anomaly and lowest en- 
ergy kink in the ARPES. Finally we found that the iso- 
tope effect on ARPES is negligible in accordance with 
experiment. These results support the claim that strong 
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